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. Abstract 

■ The electromagnetic field generated by a charged particle moving along a helical orbit 
\ inside a dielectric cylinder immersed into a homogeneous medium is investigated. Expressions 

^ ■ are derived for the electromagnetic potentials, electric and magnetic fields in the region inside 

I the cylinder. The parts corresponding to the radiation field are separated. The radiation 

■ intensity on the lowest azimuthal mode is studied. 

a\ \ 

§ ■ PACS number(s): 41.60.Ap, 41.60.Bq 

o : 

^: 

' , 1 Introduction 

Mh. 
(D . 

. The radiation from a charged particle moving along a helical orbit in vacuum has been widely 

^ \ discussed in literature (see, for instance, JIE] and references given therein). This type of electron 

motion is used in helical undulators for generating electromagnetic radiation in a narrow spectral 
interval at frequencies ranging from radio or millimeter waves to X-rays. The unique character- 
5^ \ istics, such as high intensity and high collimation, have resulted in extensive applications of this 

radiation in a wide variety of experiments and in many disciplines. These applications motivate 
the importance of investigations for various mechanisms of controlling the radiation parameters. 
From this point of view, it is of interest to consider the influence of a medium on the spectral 
and angular distributions of the radiation. 

It is well known that the presence of medium can essentially change the characteristics of the 
high energy electromagnetic processes and gives rise to new types of phenomena. Well-known 
examples are Cherenkov, transition, and diffraction radiations. In a series of papers 
the synchrotron radiation is considered from a charge rotating around a dielectric ball/cylinder 
enclosed by a homogeneous medium. It has been shown that the superposition between the 
synchrotron and Cherenkov radiations leads to interesting effects: under the Cherenkov condition 
for the material of the ball/cylinder and the particle velocity, strong narrow peaks appear in 
the radiation intensity. At these peaks the radiated energy exceeds the corresponding quantity 
in the case of a homogeneous medium by several orders of magnitude. The influence of a 
dielectric cylinder on the radiation from a longitudinal charged oscillator moving with constant 
drift velocity along the axis of the cylinder has been considered in jl2| ll.Sj. 
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In a previous paper Jl] we have investigated the radiation by a charged particle moving 
along a helical orbit inside a dielectric cylinder immersed into a homogeneous medium (for the 
radiation in a dispersive homogeneous medium see 1151). Specifically, formulae are derived for 
the electromagnetic fields and for the spectral-angular distribution of the radiation intensity in 
the exterior medium. A special case of the relativistic motion along the direction of the cylinder 
axis with non-relativistic transverse velocity is discussed in detail and various regimes for the 
undulator parameter are considered. In the present paper we consider the electromagnetic 
fields generated inside the cylinder by the charge moving along a helical orbit. The paper is 
organized as follows. In the next section, by using the Green function, the vector potential and 
electromagnetic fields are determined for the region inside the cylinder. The analytic structure 
of the Fourier components of the fields are investigated in Sec. 01 and the radiation part of the 
field is separated. A formula is derived for the radiation intensity on the lowest azimuthal mode. 
Sec. m concludes the main results of the paper. In Appendix we derive a summation formula 
over the eigenmodes of the dielectric cylinder which is used in Sec. 01 for the evaluation of the 
mode sums in the formulae of the radiation intensities. 



2 Electromagnetic potentials and fields inside a cylinder 

We consider a dielectric cylinder of radius pi and with dielectric permittivity eq and a point 
charge q moving along the helical trajectory of radius po < pi. We assume that the system is 
immersed in a homogeneous medium with permittivity £i. The velocities of the charge along the 
axis of the cylinder (drift velocity) and in the perpendicular plane we will denote by and v±, 
respectively. In a properly chosen cylindrical coordinate system (p, (j), z) with the z-axis along 
the cylinder axis, the components of the current density created by the charge are given by the 
formula 

ji = -vi6{p - po)6{(j) - uJot)Hz - v\\t), (1) 
P 

where ujq = v_l/ po is the angular velocity of the charge. This type of motion can be produced by 
a uniform constant magnetic field directed along the axis of a cylinder, by a circularly polarized 
plane wave, or by a spatially periodic transverse magnetic field of constant absolute value and 
a direction that rotates as a function of the coordinate z. In the helical undulators the last 
configuration is used. 

The vector potential of the electromagnetic field is expressed in terms of the Green function 
Gu{r,t,r' ,t') as: 

Mr,t) = I Gu{r,t,r',t')ji{r',t)dr'dt', (2) 

where the summation over / is understood. For static and cylindrically symmetric medium the 
Green function is presented in the form of the Fourier expansion 

Gii{r,t,r' ,t') = / dkzdijGii{m, kz,u>, p, p') eyip[irn{(l) — (j)') + ikz{z — z') — iuj(t — t')]. (3) 

The substitution of expressions and (jSJ into formula © gives the following result for the 
components of the vector potential: 

oo 

m=-oo •' °° 

X [v±Gi^{m,kz,uJm{kz), P, Po) + v\\Gizim,kz,u;m{kz:), p, po)] , (4) 
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where 

^^m{kz) = muJo + k^vii . (5) 

In the Lorentz gauge, using the formula for the Green function given earher in Ref. "3" and 
introducing the notations 



-^j-k;, j = o,i, 



(6) 



for the corresponding Fourier components Gu = Gii{m, kz,uJmikz), P, Po) in the region inside the 
cyhnder, /> < pi, we obtain 



IT 



l<i> 



- ^ P"' Jm+p(Ao/><)-f^m+p(AoP>) + S^''Jm+p(Aop) ,l = p,4>, (7a) 



p=±l 



Giz 



kz Jm{^OPo)Hm{^lPl) sr-^ Jm+p{^OP) rj t\ \ i 

- 1^ „a,-lwJ ^m+p{Mpi), I-P, 



PiamWm 



=±1 



m+p 



TT 

2i 



Jm{^OP<)Hm{^OP>) — Jmi^OPo)Jmi^OP) 



where Jm{x) is the Bessel function and Hm{x) = Hm\x) is the Hankel function of the first kind, 

/0< = min(p, Po)) P> = ™ax(p, po); and dp = 1, cr^ = 2. The coefficients Bm^ in these formulae 
are determined by the expressions 



W 



H 



(7b) 
(7c) 



r(p) t ^ ^m+p , iphHrn+pi^^lPl) ^ f \ ^ \ ST^ Jm+l{>^OPo) 

= - -'m+p (AoPoj 7777 \ 7— — 777-7 -nm(Aipi) 2^ 



' • 1 



m+p TrPlQ^mW^m+p 



l=±l ^^ni+l 



where 



Or 



_ iAoJ™(AoPi) y /^!!1±4^. 



m+l 



Here and below we use the notation 

= FmiXoPl) 



p dHm{XiPi) rj ^^ -.dFmiXoPi) 



dpi 



-f^m(AlPl) 



dpi 



(8) 



(9) 



(10) 



with F = J,H. Taking ei = Eq, from Q we obtain the Fourier components of the Green 
function in a homogeneous medium. In this limit = and = 0. By using formulae (jT)), 
the electromagnetic fields can be evaluated for an arbitrary motion of a charged particle with 
zero radial velocity. 

In the case of the motion along helical orbit, substituting the expressions for the components 
of the Green function into formula (j!]), we present the vector potential in the form of the Fourier 
expansion 



A{r,t)= e^-(*--oi) / dkze'''^^'-^^^'^A^i(kz,p). 



:ii) 



As the functions Ai{r,t) are real one has A'^i{kz,p) = A_mi{—kz, p)- Consequently, formula 
(jllj) can also be rewritten in the form 



Ai{r,t) = 2Re^ e^rni^-^ot) / dkze'''^'^'~^\\'^ Ami{kz, p), 

m=0 -^-^ 



(12) 
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where the prime means that the term m = should be taken with the weight 1/2. In the 
discussion below we will assume that m ^ 0. The Fourier components A^i = A„a{kz-,p) can be 
presented in the form 

^w = ^S+^S> (13) 

where corresponds to the field of the charge in a homogeneous medium with permittivity 
£0- For this part one has 

^ml = ~ Ad'y^-^ ^ p'"Jm+p{>'Op<)Hm+p{XoP>), I = pA, (14a) 

p=±l 

A^^l = ^Jm{XoP<)H„^{XoP>). (14b) 
2c 



The part due to the presence of the cylinder is determined by the formulae 

p=±i 

Agl = -'-^^Jm{X0P0)JmiX0P), (15b) 



with the coefficients 



rip)- ^^^^ Pfa) i Jm{XoPo)Hm+piXlPl)H„,iXlPl) 

cpiuimyy m^y m+p 

The electric and magnetic fields are obtained by means of standard formulae of electrodynamics. 
As is seen from formula (jllj) . analogous expressions may also be written for these fields. For the 
Fourier components of the magnetic field one has 

HUk.,p) = H^^i+H^^,, (17) 

where h'^^ corresponds to the field of the charge in a homogeneous medium with permittivity 
£0- For p > pq this part is determined by the formulae 



qXo 



p=±i 



= E ^'^^^^-(^OP), (18b) 

p=iti 



with the coefficients 



IT 



V±Jm+p{XoPo) - Vll^Jm{XoPo) 



(19) 



The corresponding expressions for p < po are obtained from (|18|) by the replacements J H 

ml 



of the Bessel and Hankel functions. The part H^j in ()17() is due to the inhomogeneity and is 
given by formulae 



= E^'''"'^- WAo/^), Z = /^,</', (20a) 



27rz'^' 
qXo 



p=±i 



= E rf^-(Aop), (20b) 

p=±i 
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where the notation 

is introduced. By making use of the Maxweh equation V x H = — iwegE/c, one can derive the 
corresponding Fourier coefficients for the electric field: 



qicXok^ 



C2 



(22a) 



where I = p,(p. As it follows from these formulae, • = 0, i.e., the corresponding Fourier 
components of the electric and magnetic fields are perpendicular to each other. 

Taking in the formulae given above ujq = for a fixed pQ, as a special case we obtain the 
fields generated by a charge moving with a constant velocity on a straight line p = po parallel 
to the cylinder axis. For this case one has 

= Af ^ - 1), = t;,! V^/c, i = 0, 1, (23) 

and v± = 0. As a result the dependence on the parameter pQ in the coefficients Dm^ is in the 
form of the Bessel function JmiXopo)- It follows from here that in the limit po ^ (the particle 
moves along the axis of the dielectric cylinder) the term with m = contributes only. The latter 
property is a simple consequence of the azimuthal symmetry of the corresponding problem. 

The formulae for the fields are simplified for the m = mode. In this case Aj = A^*^^ and for 
the function am from Q one has 

eiXoJo{Xopi)H^{XiPi) - eoXiJQ{Xopi)Ho{Xipi) 
(ei - eo)I^i 



(p) 

The expression for the coefficients takes the form 



(p) _ iTTv^Wl' t^ii Jo(AoPo) 

^0 -P o„wJ M^oPo) 



2cW( ' ' cWi 



^ttAo^// + Ho{Xipi)H'f^{Xipi) 



2k, " pi aoW^ 



(25) 



For the corresponding Fourier components of the electric field we obtain the formulae 



4/ = ^^^I^I3T--^i(AoP) JZp' ^0 > = ^^^o(AoP) \ (26) 

II " p=±i II p=±i 

with I = p,(f> and o"; is defined after formulae Q. These formulae are used in the next section 
to the derive the formula for the corresponding radiation intensity inside the cylinder. 



3 Radiation fields inside a dielectric cylinder 

In this section we consider the radiation field propagating inside a dielectric cylinder. First of 
all let us show that the part corresponding to the fields of the charge in a homogeneous medium 
with permittivity so does not contribute to the radiation field in this region. This directly follows 
from the estimate of the integral over k, in the corresponding formula on the base of the 
stationary phase method. As in the integral over k, the phase k^z has no stationary points 
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and the integrand is a function of the class C°°(i?), for large values \z\ the integral vanishes 
more rapidly than any power of From this argument it follows that the radiation field is 

determined by the singular points of the integrand in the integral over k^- In the expressions 
for the parts of the Fourier components and -E^^ , the coefficients d'^ enter in the form 
of combinations D'>m and pD^ . By using formulae (PT|) , it can be seen that for m / 
these combinations are regular at the points corresponding to the zeros of the functions 
and W^_|_^. The only poles of the parts of the Fourier components for the fields determined by 
relations (|2fl)) . (|22)) are the zeros of the function Um appearing in the denominators of Eqs. (jHJ, 
()16() . It can be seen that this function has zeros only for < 0. Note that for the corresponding 
modes in the exterior region (see Ref. [2 )> P > Pi, the Fourier coefficients are proportional 
to the MacDonald function Ky{\\i\p), v = m,mzizl, and they are exponentially damped in 
the region outside the cylinder. These modes are the eigenmodes of the dielectric cylinder and 
propagate inside the cylinder. By using the properties of the cylindrical functions, for m ^ 
the function am can also be written in the form 

""^ = (e.-eoHvi-m^u^y ^^^^ 
where we have used the notations 

K^ = \Xi\pij^ + XoPi^, u = Xo/\X,\ + \Xi\/Xo (28) 

Um = Knfeo|Ai|pi^ + eiAo/>i^Vm24^i^(eiA2 + eo|Ai|2). (29) 

Here and below it is understood Km = Km{\Xi\pi), Jm = </m(AoPi) if the argument of the 
function is omitted, and the prime means the differentiation with respect to the argument of 
the function. Now the equation for the eigenmodes is written in the standard form (see, for 
instance, [T6| ) 

Um = 0. (30) 

Unlike to the case of the waveguide with perfectly conducting walls, here the eigenmodes with 
m / are not decomposed into independent TE and TM parts. This decomposition takes place 
only for the m = mode and this case will be discussed below separately. 

We denote by Aopi = Xm,s, s = 1,2, . . ., the solutions to equation pOf) with respect to Aopi 
for the modes with m ^ 0. The corresponding modes kz = km}s are related to these solutions 
by the formula 



U±) 



c 1-/3; 



l + &ls(/3o'|| 



cA™,.,, ^^^^ 



These modes are real under the condition 6^ ^1 — /9q|| j ^ 1. For /3o|| < 1, the condition for 
km}s to be real defines the maximum value for s, which we will denote by Sfyi'. 



If the Cherenkov condition /3o|| > 1 is satisfied, the upper limit Sm for s is determined by the 
dispersion law for the dielectric permittivity via the condition eo{iOm) > c^/^u- The function 
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^m{kz) is also expressed in terms of 



1 



Pi. 



1 ± 



1 + 62 



(33) 



In order to obtain unambiguous result for the fields we should specify the integration contour 
in the complex plane kz in Eq. For this we note that in physically real situations the 

dielectric permittivity eo is a complex quantity: Eq = + *£o- Assuming that Eq is small, this 
induces an imaginary part for kz given by formula 



Im k^ 



l±/3o||,/l + fe^,,(/32 



1^ 



e^(l-/52|)2^1 + 6^^^(/3, 



(34) 



where the coefficient and the argument of eg(u;m) ai^e evaluated with the real part of dielectric 
permittivity. Note that one has eo('^m) ^ for ^ 0. In the case /3o|| < 1 from we 
have ujm > and, hence, in accordance with (|34|) . ImA;^ ^ for the modes kml. Deforming the 
integration contour we see that for /3o|| < 1 in the integral over kz in H12|) . the contour avoids the 

poles kmji from above and the poles kml from below. For the case /3o|| > 1 we have LOm < for 

the modes with kz = km^s and ujm > for the modes with kz = fcm,l and, hence, in accordance 
with formula 1)341) for both types of modes Imkz < 0. As a result, in this case in the integral 
over kz all poles km} should be avoided from above. 

Now let us consider the electromagnetic fields inside the cylinder for large distances from the 
charge. In this case the main contribution into the fields comes from the poles of the integrand. 
Having specified the integration contour now we can evaluate the radiation parts of the fields 
inside the cylinder. Closing the integration contour by the large semicircle, under the condition 
/9o|| < 1 one finds 



cr47rRe 



oo 
. m=l 



im{(f>—uJot) 



V Res e^'^'^'-^ll'^Fmiikz^p) 
^ I, _i.(<^) 



s=l 



(35) 



where a = +(— ) for z — v^^t > {z — v^^t < 0). This expression describes waves propagating along 
the positive direction of the axis z for a = + and for a = — , bm,s < and waves propagating 

along the negative direction to the axis z for o" = — , 1 < bm,s < 

condition /3o|| > 1 for the radiation fields one finds 



'1 



/32||. Under the 



F['^''\r,t) = -^F^''\r,t)eiv^^t 



(36) 



a=± 



where 9{x) is the Heaviside unit step function. As we could expect, in this case the radiation 
field is behind the charge. Formula ()36|) describes the waves propagating along the positive 
direction of the axis z for a = + and for a = — , bm,s > 1, and waves propagating along the 
negative direction of the axis z for a = — , bm,s < 1- Note that for bm,s > 1 there are no waves 
propagating along the negative direction of the axis z. 

By evaluating the residue in ()35)) . for the functions F^'^\Y,t) one finds 



Fl''\v,t) 



47rRe 



m=l 



E 

s=l 



im{4>—uiot) \ ikziz~vat) -^ml 



dam/dkz 



(37) 
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where the formulae for the coefficients pjf') [kz, p), I = p, (p, z are obtained from the corresponding 



expressions for by the replacement Cm^ — > Cm"^^ for the components of the vector potential 
and by the replacement Dm^ Cm'^^ for the electric and magnetic fields, with 



K K 



A hi^) Jm{^m,sPOl) ^(o■) Jm+l{^m,sPOl) 



(38) 



where we have introduced the notation 
and 

m = J (l - k'^s'pl - JiA^,, ,01 = ^. (40) 
y V ^0/ £0 Pi 

Note that we have the following relations 

2 

\ X^*^) T/ 



Wl^lp = {pVm - mu) JmKra, wj^^ = ^ ^"'''[^■^\ ^ [V^ - mW) J^K^. (41) 



Unfortunately, for the case m 7^ we have no summation formula over the eigenmodes like 
to that for the m = mode (see below) and for the further evaluation of the radiation fields 
numerical methods have to be used. 

Now we turn to the modes with m = 0. In this case there are two types of eigenmodes. The 
first ones are the zeros of the function ao corresponding to the zeros of the nominator in (|24|) . 
The equation for these eigenmodes is written in the form 

^o|Ai|^ + eiAo|^ = 0, A, = Af , (42) 

where A^^"* is defined by formula (|23|) . For the presence of these modes we have the conditions 
(3i\\ < 1 < /3o||. In particular we should have ei < eo- From formula (|25|) it follows that the 

combination "^ppDm^ has no poles at these eigenmodes and for the corresponding radiation 

field one has H^J = H^J = and E^J =0. As a result, these waves are TM waves or the waves 

of the E-type. The values for kz defined by equation (|^^ we will denote by ^k^] , s = 1,2, . . .. 
The second type of m = eigenmodes corresponds to the zeros of the function Wf . The 
corresponding equation can be written in the form 

Ao^ + |Ai|^ = 0, A,, = Af. (43) 
-Ko Jo 

These zeros are poles for the combination pDm^ , whereas the combination Dm^ is regular 

at these zeros. For the corresponding radiation fields one has H^]^ = and E^^J = E^^J = and 
these waves are TE waves or the waves of the M-type. The corresponding eigenvalues for kz we 
will denote by zL/cq^^ s = 1,2, . . .. Note that both types of the eigenmodes are also obtained 
from H3U|) taking m = 0. 

As in the case of m 7^ 0, to define the fields from formula (|12|) we should specify the contour 
for the integration in the complex plane kz- In the way similar to that used for the modes with 
m 7^ 0, it can be seen that the contour avoids the poles itfcg J from above. The corresponding 
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radiation fields are found by evaluating the /c^-integrals with the help of residue theorem. For 
z — v\\t > we close the integration contour in the upper half-plane of the complex variable 
kz and the radiation fields are zero. We could expect this result as /?o|| > 1 and the radiation 
field is behind the charge. For the points in the region z — ■uyt < we close the integration 

(V) 

contour in the lower half-plane and the integral is equal to the residues at the poles iL/cq J and 



(M) 



g multiplied by —2iri. Here we will present the formulae for the nonzero components of 
the electric field: 



E 



(E) 

p,m=0 



Aqy. 



^1^^^ Jo(/^oAEVi(pAi5) sin 



^0,s ) 



E 



(E) 

z,m=0 



4g 



1) Mpo>'oJ)MpKJ)cos 



V\\t) 



£1 — eo 



E 



(M) 
(f),m=0 



iqv± ^ Ji{Po\^s )JiiP\ 



eoJl?{piy^J)+ei{f3',^^ - l)Ji{piX'o7) 



ph 



(ei - eoKHpi^Z^) 
where we have introduced the notation 



cos 



viit) 



A 



(F) 
0,s 



1, F = E,M. 



(44) 



(45) 



In these formulae the summands with a given s describe the radiation field with the frequency 



(F) 
0,s 



viik 



(F) 
0,s • 



In Appendix it is shown that Xq^Pi ^ 1 and, hence, for the corresponding 



frequencies one has the estimate u>l^^ > v\\/{piyj P^^^ — 1). In formulae H44() the upper limit for 
the summation over s, which we will denote by Sm, is determined by the dispersion law of the 
dielectric permittivity e = e{uj) through the condition e{ujl^]) > c^/fy ■ 

Having the electric field we can evaluate the radiation intensity on the mode m = inside 
the dielectric waveguide by using the formula 



Im=0 

For this quantity one finds 



{U^<f>,m=o + jzEz,m=o) pdpdifdz. 



r{E) 
'm=0 



r{M) 
'm=0' 



lm=0 

where for the contributions of the TM and TM waves we have 



(E) 
m=0 



AM) 
'm=0 



Pl^\ 



E 

s 

E 



1 



'^oiPO>%s) 



Eq — Si 



)^ 

0,s 



{E)^ 



JHpoX^^ 



(£o-^i)J^2(piAi'fV 



(46) 
(47) 

(48a) 
(48b) 



As the modes kj^^ are not explicitly known as functions on s, formulae (|48a() . H48b|) are not 
convenient for the evaluation of the corresponding radiation intensities. More convenient form 
may be obtained by making use of the summation formula (|58|) derived in Appendix, taking 



1 - P'n 



(49) 
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In formula (|58j) we choose a = Eq/ei, f{z) = zJq{zpq/ pi) for the waves of the E-type and 
a = 1, f{z) = zJf{zpo/pi) for the waves of the M-type. For both types of modes one has 
f{—ix) = —f{ix) and we can use the version of the summation formula given by (|6U|) . In 
the intermediate step of the calculations, it is technically simpler instead of considering the 
dispersion of the dielectric permittivity to assume that in formulae H48() a cutoff function V^^CAg^j) 
is introduced with p being the cutoff parameter and -00 = 1 [for example, ip^{x) = exp(— /xx)], 
which will be removed after the summation. In this way, after the application of the summation 
formula we find the following results 



(E) 



m=0 



r(M) 
'm=0 



g^f II 



2 2 
g f]_W|| 



dx—'il)f,{x)jQ{pox) + 



dx 



£0 

x'4)^{x) 



tt'^pI 



dx 



Ilixpo/pi 



^iXQeo/eA^^^) 



1 



jlipQx) 



^^p\ 



dx 



Ilixpo/pi) 



(/3oii - '^)xgi{'n,x) 



(50a) 



(50b) 



where the function gai^l^x) is defined by formula ()61() . In the first terms in the square brackets 
replacing the integration variable by the frequency uj = v\\x / ^ (i^^^ — 1 and introducing the 
physical cutoff through the condition /Jgy > 1 instead of the cutoff function, we see that these 
terms coincide with the radiation intensities on the harmonic m = for the waves of the E- and 
M-type in the homogeneous medium with dielectric permittivity eq (see Ref. PH|)- The second 
terms in the square brackets are induced by the inhomogeneity of the medium in the problem 
under consideration. Note that, unlike to the terms corresponding to the homogeneous medium, 
for pq < pi the latter are finite also in the case when the dispersion is absent: for large values 
of x the integrands decay as exp[— 2x(l — po/pi)] (for this reason we have removed the cutoff 
function from these terms). In particular, from the last observation it follows that under the 
condition (1 — po/ pi) ^ v\\/uJd-, where uJd is the characteristic frequency for the dispersion of 
the dielectric permittivity, the influence of the dispersion on the inhomogeneity induced terms 
can be neglected. Note that in a homogeneous medium the corresponding radiation propagates 
under the Cherenkov angle 9c = arccos(l//?o||) and has a continuous spectrum, whereas the 
radiation described by H50|) propagates inside the dielectric cylinder and has a discrete spectrum 
with frequencies i^^} ■ As the function gaiv^^) is always nonnegative, from formulae (|5()j) we 
conclude that the presence of the cylinder amplifies the m = part of the radiation for the 
waves of the E-type and suppresses the the radiation for the waves of the M-type. In the helical 
undulators one has v± <C iJy and the contribution of the TM waves dominates. 

In addition to the total intensity, it is of interest to have the distribution of the intensity 
over the frequencies. To this aim we define the partial intensities Z^=o s fo'^ ^ given s related to 



the quantity I'n by the formula / 



(F) 
m=0 



Z-^s m=0,s ■ 



In figure n we have presented the results 



of the numerical evaluation for the number of the quanta radiated per unit time, N 



(F) 



r(E) 
m=0,s 



(V) 

/hujQ J , on a harmonic with a given s, as a function of s for the electron with the energy 
2 MeV moving inside a cylinder with dielectric permittivity eq = 3. It is assumed that for 
the exterior region ei = 1. The angle between the electron velocity and the cylinder axis 
(tan^o = t^±/t^\\) is taken 0.1 rad, and po/pi = 0.5. The values of the parameter PiAqj, for the 
first eight harmonics are presented in the table. The corresponding frequencies are related to 
this parameter by the formula oJq J 



m=0,s 



1. 



s 


1 


2 


3 


4 


5 


6 


7 


8 




2.728 


5.931 


9.102 


12.262 


15.415 


18.565 


21.713 


24.860 




2.507 


5.653 


8.801 


11.947 


15.092 


18.236 


21.380 


24.523 
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r i r I , 



t t t 



Figure 1: The number of the radiated quanta on the mode m = as a function of the harmonic 
number s for the waves of the E- (left panel) and M-type (right panel) (for the values of the 
parameters see the text). 

4 Conclusion 

We have considered the electromagnetic field generated by a charged particle moving along a 
helical orbit inside a dielectric cylinder immersed into a homogeneous medium. The fields and 
the spectral-angular distribution of the radiation intensity in the exterior medium have been 
investigated in our previous paper particular, the conditions were specified under which 

strong narrow peaks appear in the angular distribution for the number of radiated quanta. In 
the present paper we study the fields inside the cylinder. By using the corresponding formulae 
for the components of the Green function, we have derived expressions for the electromagnetic 
potentials and fields. These expressions are presented in the form of sums of two terms. The first 
ones correspond to the fields generated by the charge in a homogeneous medium with the same 
dielectric permittivity as that for the material of the cylinder. The second terms are induced by 
the difference of the dielectric permittivities for the exterior and interior media, and are given by 
formulae (|15|) , H20|) , 1)22(1 . We have extracted the parts of the fields which are responsible for the 
radiation. For the radiation propagating inside the cylinder the projection of the wave vector 
on the cylinder axis takes discrete set of values which are solutions to the eigenmode equation 
()3fl|) . The expressions for the electromagnetic fields corresponding to the radiation propagating 
inside the cylinder are derived. As an application we have considered the radiation intensity on 
the mode m = 0. In this case the radiation field is separated into purely TE and TM modes 
and we have found the intensities for both types of these modes determined by formulae H48|). 
These formulae presents the intensities as sums over the modes of the dielectric waveguide. 
However, as the corresponding eigenmodes are not explicitly known, this form is not convenient 
for the numerical evaluation of the radiation intensity. In Appendix, by using the generalized 
Abel-Plana formula, we derive a summation formula for the series over the normal modes of the 
dielectric cylinder which allows to extract from the radiation intensity the parts corresponding to 
the radiation in a homogeneous medium and to present the inhomogeneity-induced part in terms 
of exponentially converging integrals. Unlike to the parts corresponding to the homogeneous 
medium, the terms induced by the presence of the cylinder are finite also in the case when the 
dispersion is absent. We have specified the condition under which the influence of dispersion on 
the inhomogeneity-induced effects can be neglected. 
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A Summation formula over the eigenmodes 

In this section we derive a summation formula for the series over zeros of the function 

C^{n,z) = Va{Mz),K^{r^z)], (51) 

where and in what follows for given functions F{z) and G{z) we use the notation 

V{F{z),G{r]z)] = F{z)G'{r]z) + arjF' {z)G{riz), (52) 

with a ^ 1 and r] being real constants. For this we will use the generalized Abel-Plana formula 
from (for applications of the generalized Abel-Plana formula in quantum field theory with 
boundaries see JHl)- In this formula we substitute 

ya{Yo{z),Ko{vz)} 



g{z) = if{z)- 



(53) 



Ca(ry,z) 

For the combinations of the functions entering in the generalized Abel-Plana formula one has 

y^{H^^\z),Ko{vz)] 



f{z)-{-lfg{z)=f{z)- 



(54) 



C«(r?,z) 

where Hi^\z), k = 1,2, are the Hankel functions. Let us denote positive zeros of the function 
Caif], z) by kg, s = 1,2, . . ., assuming that these zeros are arranged in the ascending order. Note 
that, for z » 1 we have Ca{r],z) « K'Q{rjz) — arjzKQ{r]z)/2 < and, hence, kg > 1. By using 
the asymptotic formulae for the cylindrical functions for large values of the argument, it can be 
seen that for large values s one has fcg ~ — arctan(l/ar/) + tt/A + tts. For the derivative of the 
function Ca{r],z) at the zeros kg one obtains 

K^ivz) 



Mz) 



[a{l + ar]'^)Jo^{z) + {a 



l)Jl{z) 



ks 



(55) 



In particular, it follows from here that the zeros are simple. Assuming that the function f{z) is 
analytic in the right half-plane, for the residue term in the generalized Abel-Plana formula one 
finds 



liesg{z) 

where we have introduced the notation 

Ta{z) 



-T^ikg) f {kg) , 



2a/z 



(56) 



(57) 



a{l + ar]2)4^(z) + (a - 1) J2(z) ' 

Substituting the expressions for the separate terms into the generalized Abel-Plana formula we 
obtain the following result 



lim 



■^0 nxo 1 /" 

TTo,{kg)f{kg)- / dxfix) = — / 



dx 



K.{i^o(x),fff (r?x)} 
VMx),H^i\ilx)] 
(1)/ 



+fi-ix 



Va{Ko{x),H;,''{vx)} 
Va{Io{x),Hl,'\r,x)} \ 



(58) 
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where sq is defined by the relation k^g < xq < kg^^i. This formula is valid for functions f{z) 
obeying the condition 

\f{z)\ < e(x)e''l^l, z = x + iy, \z\ oo, (59) 

where c < 2 and e{x) — > for x oo. Formula H58|) is further simplified for functions satisfying 
the additional condition f{—ix) = —f{ix): 



where 



lim 



.3 = 1 



dxf{x) 



Aia , f(ix) 



Jo x'^ga{v,x)' 



(60) 



Il{x) [Jf (ryx) + Y^i'qx)] + a^llix) [Jl{r,x) + Y^ir^x)] 
-2ar]Io{x)Ii{x) [Jo{r]x)Ji{r]x) + Yo{r]x)Yi{r]x)] . 



Note that we have denoted ga{'q,x) = \V{Iq{z), H''^\riz)^\'^ and this function is always nonneg- 
ative. 



(61) 
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